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We propose and analyze a scheme to create dark solitons 
and vortices in Bose-Einstein condensates. This is achieved 
starting from a condensate in the internal state \a) and trans- 
ferring the atoms to the internal state |6) via a Raman tran- 
sition induced by laser light. By scanning adiabatically the 
Raman detuning, dark solitons and vortices are created. 
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Recently Bose-Einstein condensation has been demon- 
strated in dilute atomic gases 0. This state of matter 
resembles other states found in the fields of superfluidity, 
superconductivity and nonlinear optics. It is thus natu- 
ral to expect that some phenomena that appear in those 
fields, such as solitons and vortices, can be observed with 
dilute atomic gases. In fact, the Gross-Pitaevskii Equa- 
tion (GPE) Q which describes the wavefunction of the 
macroscopically occupied state in a trapping potential al- 
lows for stationary solutions that represent dark solitons 
(see below) and vortices (see also ||). Here we propose a 
scheme to generate these solutions in a controlled way us- 
ing an approach based on "engineering the macroscopic 
wavefunction" : once the condensate has formed, we use 
a coherent Raman process to create solitons and vortices. 

Our idea is to couple the internal state \a) where the 
condensate is formed with another internal state \b) us- 
ing a Raman transition (Fig. 1). The laser parameters 
are chosen such that the state after the transfer is an 
eigenstate of the GPE corresponding to solitons or vor- 
tices. The complete transfer is achieved by an adiabatic 
change of the laser frequency. In the present case, where 
the whole process is described by non-linear equations, 
the familiar idea of adiabatic transfer along eigenstates 
has to be taken with caution since the Hamiltonian de- 
scribes the atomic interactions by a mean field depending 
on the shape of the wavefunctions. We will analyze two 
cases: in ID we will study the creation of dark solitons 
and in 2D the creation of vortices. 

In the Hartree-Fock approximation the state <!> of a 
condensate of N bosons confined in a potential V(r) is 
described by the time-independent GPE 
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The mean field interactions are characterized by a cou- 
pling constant g — AirTi 2 a s /m, where a s > is the s-wave 
scattering length. In the following analysis we will con- 
centrate on the Thomas-Fermi limit H] , since most of the 
experiments operate in this regime m. In this limit, the 
mean interaction energy is much larger than the mean 
kinetic energy which can be neglected when calculating 
the ground state solution of the GPE (fil) : 



$ Eo (r) = {[E -V(r)]/(Ng)} 1 /\ 
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for f such that V(r) < Eq: Eq is determined from 
J d 3 r\^E ('f f )\ 2 = 1 which reflects particle conservation. 
Here we are interested in other stationary solutions of 
(0) <j>(r) with energy larger than Eq. Far from the trap 
center r = we can still neglect the kinetic energy which 
suggests the ansatz 



$(f) = 0(r)$ s (r) 
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that is a a product of an envelope function (g) and a func- 
tion (j) with the condition |0(r)| ~ 1 far away from the 
origin. Now E is determined from J d 3 r\4>(r)$E(f)\ 2 = 1. 
Substituting (||) in the GPE and neglecting the deriva- 
tives of $_E(r) we obtain a nonlinear equation for <j)(r), 
which near r — 0, where we can neglect the variation of 
the trapping potential, reads 
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This equation has the same form as the familiar GPE 
for the homogeneous case. Together with the boundary 
conditions |<^(r)| ~ 1 for \r\ — > oo it gives rise to dark 
solitons and vortices 0. Our goal is to design boson- 
laser interactions which will generate these solutions. 

We assume that the bosons have two internal levels \a) 
and | b) as in Rb ||. The particles interact with a laser 
beam that connects these two levels. The evolution of 
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the mean field spinor $ = (<!> a , <!>;,) obeys the following 
non-linear equation 
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where 
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+ V(f) + Ng[\<P a (r,t)\ 2 + |$ 6 (r,i)l1, (6) 



describes the linear evolution plus atom-atom interac- 
tions and 
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the interaction with the laser. In (g) we have assumed 
that the interaction between levels i and j (i,j = a,b) 
can be described by a pseudopotential gijS(f) and that 
9ij = 9 > Oj which indeed is a good approximation 
for Rubidium 0. For V(f^ we choose an anisotropic 
harmonic potential with frequencies u> x . V:Z which we as- 
sume identical for both internal levels. In ([?]), A(r) = 
f2 a (r)f2b(r)/(4A) where A is the detuning from the in- 
termediate level |r), and tt a ,b are the Rabi frequencies 
corresponding to the couplings between \a) and \b) to 
|r), respectively (see Fig. 1). Their specific form de- 
pends on the laser configuration. The Raman two- 
photon detuning is denoted by 8. The conservation of 
the number of particles gives the normalization condi- 
tion Jd 3 f [|$ a (r)| 2 + IM^l 2 ] = *l + A = 1 with nl b 
the populations in levels a, b, respectively. 

As the initial state we take &e [see Eq.(||)] which cor- 
responds to the state of the condensate formed in the in- 
ternal level | a). We will design A(r) and 8 such that the 
atoms are transferred to |6) with a wavefunction which 
corresponds to dark solitons or vortices. In absence of 
interactions (g = 0) the problem reduces to the one of 
a single trapped particle. In that case, one can sim- 
ply use a resonant 5 = laser pulse of a well defined 
area to carry out the population transfer M. In pres- 
ence of interactions, this method will not work: as soon 
as particles are transferred to a different state and the 
shape of the wavefunction changes, the interaction en- 
ergy changes [Fig. l(a,b)]. Therefore, an initially res- 
onant Raman laser pulse at the beginning of the pulse 
soon becomes off-resonant, and the transfer process will 
stop. We circumvent this problem by using adiabatic 
passage. The idea is to start from a negative Raman de- 
tuning so that the atoms do not feci the laser [Fig. 1(a)]. 
Then, the Raman detuning is changed adiabatically to 
sufficiently large positive values [Fig. 1(b)]. As soon as 
the laser frequency approaches the Raman resonance, the 
atoms will start flowing to the state \b). The fact that 
the interaction energy changes will effectively change the 
value of (5 at which this resonance occurs. This will not 
affect the overall process provided the final value of 5 is 



large enough so that at the end the atoms do not feel the 
off-resonant laser anymore. The reason is that adiabatic 
transfer only depends on the initial and final values of 
the adiabatic parameter (detuning). 

In order to describe analytically the adiabatic process 
we look for stationary solutions $a(r) of (||) for a given 
value of 8. The idea is to change 8 adiabatically so that 
the state of the system changes according to in accor- 
dance to the adiabatic theorem. We thus have to impose 
that for the initial and the final Raman detunings 8q 
and Sf the spinor $5 corresponds to the initial state and 
the desired state, respectively. That is: (i) for 8 — ► So, 

$s - ($£ >°); (ii) for 6 ~> s fi *<5 = (0,$fc) where $ h is 
another stationary solution of Eq. (|l|). In order to find 
the appropriate stationary solutions of (|5|) we have to 
solve first the set of nonlinear equations 



(8) 



for a fixed (real) value of n a (with 7Tf, = |1 — it 2 ] 1 / 2 due to 
normalization) and satisfying the conditions (i,ii). Here 
^o°6 describe the wavefunctions in absence of laser cou- 
pling for given populations 7T 0i b. The presence of the laser 
will lead to a dressing of these levels in analogy with the 
well-known picture known from non-interacting atoms in 
presence of laser light M . Once these functions are found, 
we replace in (||) by (||), that is we restrict the evo- 
lution to the subspace defined by these two functions; by 
multiplying the first equation by $ a (f)* and the second 
by $b(r)* and integrating we obtain 
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Equations fl9) define two generalized dressed states of our 
system which take fully into account atomic interactions. 
From these equations we can determine the values of 8 
and E corresponding to $„ a 6 . In summary, the problem is 
reduced to solve the coupled eigenvalue equations (||) for 
a given value of 7r a , such that varying continuously this 
parameter we go from Q a {r) = ^b (' 5 ) to the desired 
state $fc(r). The energy separation A a: fc(7r a ) = e^ a — ej a 
between the bare wavefunctions <&Z°b gives the resonance 
condition for the Raman detuning, 8 — > A a j,(7r a ), for 
fixed level populations 7r ai 6. The induced width is given 
by A(7r a ) for ir 2 . ~ 1/2 which gives the avoided crossing 
and therefore the time scale for adiabaticity. 

We illustrate this procedure now for the ID case. This 
corresponds to the limit in which ui x . y 3> oj z so that the 
dynamics along the x and y direction is frozen. Our goal 
is to create a dark soliton starting from <I>£ (z). We are 
interested in dark solitons with a zero at the trap center. 
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This requires that the laser interaction changes the par- 
ity of the wavefunction when the atoms are transferred 
from a to b. To this aim we choose the simplest laser 
configuration, so that A(z) = Aosin(fcz), i.e. a stand- 
ing wave. In order to achieve an efficient coupling we 
take k < 1/zq, where zq is the size of the <J>£ ||; note 
that the effective A defined in ( pp| ) will be very small if 
kz 3> 1. In this case the avoided crossing of the dressed 
energy levels will be of the order A , which sets the time 
scale for the adiabaticity. On the other hand, Ao has 
to be smaller than the typical energy separations \A a b\ 
so that the Stark shifts do not mix these wavefunctions 
with others of higher energies e. The initial value of the 
Raman detuning must be Sq <C A a ^(ir a = 1), whereas 
the final value must fulfill Sf ^> Aa^^a = 0). In Fig. 2 
we have plotted numerical results of the solutions of the 
time-dependent GPE §). Figure 2(a) shows the spatial 
distribution P a ,b( z ) = l^a.&WI 2 corresponding to states 
| a) (solid line) and \b) (dashed line). As the transfer pro- 
gresses, we see that the wavefunction of the atoms in \a) 
narrows and the one of the atoms in \b) develops a hole 
in the center. When most of the atoms are in b, the cor- 
responding wavefunction contains a dark soliton. This 
manifests itself in the effective (trap plus mean field) po- 
tential; it is initially flat, and later it develops a narrow 
dip as a consequence of the dark soliton [see Fig. 1(b)]. 
The part of the atoms still in \a) become trapped in a 
bound state of this dip, which becomes deeper as we move 
more atoms in the excited state. In Fig. 2(b) we have 
plotted the fraction of atoms P a f, = & of levels a and 
b. As this figure shows, the transfer efficiency is essen- 
tially 100 percent. For the analytical understanding of 
these results, we proceed as explained above in terms of 
generalized dressed states. First, in the limit n a — > 1 
one can estimate the value of A a b using a square well 
of length equal to the size of the Thomas Fermi solution 
$£; . In the opposite limit, n a — > we can calculate A a b 
in the same way as in the context of Eq. (0) . We write 
$(z) = [4> a (z),<f) b (z)}<£> E (z), with $ B (z) defined in (g). 
Near the trap center we obtain 
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<f>a,b( z ) = £ a ,b4>a,b{z), 
(11) 



with the boundary conditions |</> a (z)| — > and |0b(z)| — > 
1 as \z\ — > co. We find that <t>a{z) = Asech(az) and 
<fib{z) = tanh(az) solve Eq. (|TT|). With the normalization 
condition we can find the values of A, a and t a ,b for a 
given value of 7r a . This solution is in perfect agreement 
with our ID numerical results. We performed a full 3D 
integration of the GPE in order to make sure that the 
ID effect is not effected by the presence of the transverse 
degrees of freedom |n| . 

A two-soliton solution is obtained by starting from the 
ground state in a, and coupling with a laser configuration 



which preserves the parity, X(z) = Xocos(kz). In order 
not to couple to the ground state in b, the initial detuning 
has to be #o > |Ao|. We then increase the detuning adia- 
batically to a sufficiently large value. In Fig. 3 we show 
plots of numerical solutions. At the end of the process all 
the particles are in the state \b) with a wavefunction that 
includes two dark solitons. Again an analytical Ansatz 
is possible: for |-7T£, | — * 1 we set $6(2) = 4>b(z)§E(z) with 
4>b{z) — tanh[a(z — a)] tanh[a(z + a)] (a is a free param- 
eter) which reproduces the numerical results very well. 

A 2D situation arises in the limit lu z 3> = uj x = 
Lu y . In this limit we are interested in creating vortex 
solutions of the form <&(/?, ip) = f(p)e lv where p and <p 
are cylindrical coordinates and f(p) is a function with a 
zero at p = 0. In order to provide the required angular 
momentum to the atoms that are transferred, we choose 
a laser configuration such that X(x,y) = Xo[sin(kLx) + 
ism(k L y)] ~ X k L pe iV for k L p < 1 @. The density 
distribution |<I>(x,?/)| 2 after an adiabatic switch of the 
detuning is plotted in Fig. 4. The insert shows that all the 
population is transferred to the vortex state. Analytical 
approximations can be obtained with the ansatz <p a {f) = 
Asech(ap) and 0t(r) = tanh(ap)e lv . 

A simple way of observing the shape of the density n(f) 
is by opening of the trap. In a way similar to Ref. [l2| | we 
can show that the density at later times is related to the 
density at t — by n(r,t) — n[r/^{t),Q}/^ 2d {t) where 
the scaling factors obey A = oj 2 /X d+1 (d is the dimen- 
sion) and initial conditions A = 1,A = 0. This leads to 
an asymptotic behavior 7(i) — > \[2uj z t for the ID dark 
soliton and j(t) — > uj±t for the vortex solution. This self- 
similar expansion without a change in shape is typical for 
solitonic behavior. On the other hand, another impor- 
tant issue to address is the stability of vortices and dark 
solitons. In some recent works dealing with excitations 
of vortex states H it is shown that quasiparticles states 
localized near the center of the trap will be preferentially 
occupied by collisions thereby destabilizing the vortex. 
However, using the analogy between our proposal and 
the excitations of condensates via time-dependent trap- 
ping potentials [OL we expect the destabilization time 
to be much longer than that required for the creation of 
vortices and solitons. The reason is that the creation of 
vortices and dark solitons involve regular (non-chaotic) 
solutions of the GPE, which give rise to destabilizations 
that grow only polynomially with time (instead of expo- 
nentially) ||. 

We have demonstrated that one can engineer the 
macroscopic wavefunction of Bose-Einstein condensed 
sample by coupling the internal atomic levels with a laser. 
The method is based on adiabatic transfer of population 
along generalized dressed states which include the non- 
linear atom-atom interactions. We find the technique to 
be very robust against uncertainties in the parameters of 
the problem. Furthermore, we have developed analytical 
approximations to describe this process, and show with 
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explicit examples how to generate dark solitons and vor- 
tices. Our numerical results confirm our predictions and 
demonstrate the stability of these solutions. We expect 
that this method can be applied to current or planned 
experiments. 

We thank K. Burnett for discussions. This work was 
supported in part by the TMR network ERB-FMRX- 
CT96-0087, and by the Austrian Science Foundation. 



FIG. 3. Generation of two dark solitons. The detuning is 
varied linearly with time from 8 — 0.25a; to S = 5o>. Other 
parameters: Ao = 0.15a;, kh = 0.5/ao and Ng = 207kjao, 
where ag = %/{muj). 



FIG. 4. Generation of a vortex: position distribu- 
tion of the final state. The detuning is varied linearly 
from 8 = — 0.6cjx to 5 = 5lo±. The parameters are 
Ng = 500hu>±a 2 L , A = 0.157ia;_i_ and kh — 0.5a±. The in- 
set shows the evolution of the populations in levels |o) and |6) 
(solid and dashed lines, respectively). 
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FIG. 1. Schematic representation of the process: (a) ini- 
tial state; (b) final state. 



FIG. 2. Generation of a dark soliton. The detuning is var- 
ied linearly with time from 8 — —1.5a; to 8 — 6.5uj. Other 
parameters: Ao = 0.15a;, kh — 0.5/ao and Ng — dOTuoao, 
where a§ = %/{muj). (a) snap shots of the position distri- 
butions of the wavefunctions corresponding to atoms in level 
| a) (solid line) and |6) (dashed line) for different times; (b) 
Populations of these levels as a function of time. 
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